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ABSTRACT 

Experimental and computational studies have been 
performed to investigate the detailed distribution of convective 
heat transfer coefficients on the first-stage blade tip surface for a 
geometry typical of large power generation turbines(>100MW). 
In a previous work the numerical heat transfer results for a 
sharp edge blade tip and a radiused blade tip were presented. 
More recently several other tip treatments have been considered 
for which the tip heat transfer has been measured and 
documented. This paper is concerned with the numerical 
prediction of the tip surface heat transfer for radiused blade tip 
equipped with mean-camberline strip (or “squealer"’ as it is 
often called). The heat transfer results are compared with the 
experimental results and discussed. The effectiveness of the 
mean-camberline strip in reducing the tip leakage and the tip 
heat transfer as compared to a radiused edge tip and sharp edge 
tip was studied. The calculations show that the sharp edge tip 
works best (among the cases considered) in reducing the tip 
leakage flow and the tip heat transfer. 

NOMENCLATURE 

Cp constant pressure specific heat 
h heat transfer coefficient 

Pr Prandtl number 

R gas constant 

Re Reynolds number 

T static temperature/ 

Tu turbulence intensity 

V magnitude of the velocity/(/? Tq) JTTp 

dimensionless distance from a wall, = — 

Y specific heat ratio 

X shear stress 


Total pressure loss. 




Subscripts 

t total conditions 

w wail value 

0 inlet condition 

1 exit conditions 


INTRODUCTION 

Blade tips are susceptible to burnout and oxidation due to 
high thermal loading caused by the tip leakage flow. Efficient 
internal or film cooling schemes are necessary to protect turbine 
blades against damage. The design of such schemes requires 
detailed knowledge of heating patterns on and near the tips which 
could be gained by experimental and predictive methods. It was 
shown by Mayle and Metzger [ 1 ] that the heat transfer at a blade 
tip (and by extension total pressure loss) is not very much affected 
by rotation. This has spurred experimental work on measuring tip 
heat transfer for various blade tip treatments on stationary 
cascades and has provided data for numerical simulations. In a 
previous numerical study [2] the predicted tip heat transfer for a 
sharp edge tip and a radiused tip (blade tip with a rounded edge 
around the perimeter of the blade tip) was compared with the 
measured experimental data of Bunker et al. [3]. The data was, at 
the time, the only set of liquid crystal measurements of tip heat 
transfer. The results of simulations nicely agreed with the 
experimental data and showed that the tip heat transfer for 
raffiused edge would be larger than that for a sharp edge tip due to 
the relatively higher velocity on the radiused edge tip than on the 
sharp edge tip. More recently. Bunker and Bailey [4] have made 
further measurements for various tip treatments including 
circumferential rub strip, 45 -degree angled run strips, perimeter 
squealer rim and mean-camberline strip. 
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The perimeter squealer tip case has already been computed 
by Ameri et al. [5] though for a different blade configuration. In 
this paper the mean-camberline strip (squealer) case has been 
computed. The computed tip heat transfer has been compared 
with experimental measurements. As the experiment was mainly 
concerned with the measurement of the rate of heat transfer, the 
numerical simulation attempts to construct the flow structure in 
order to help explain the phenomena observed in the 
experiments. 

In the ensuing section a brief description of the experiment 
simulated in this work is provided. The numerical method 
employed and the grid topology used to model the 
computational domain are described subsequently. The results 
of the simulations showing the flow characteristics and 
computed tip leakage rate as well as the tip heat transfer are next 
presented and are followed by a summary and presentation of 
the conclusions of the work. 

THE EXPERIMENTAL SETUP 

The experimental setup and conditions are similar to those 
described in [3]. In that paper the pressure distribution and tip 
heat transfer for both a sharp edge tip and radiused edge tip of 
the present blade geometry were presented. The heat transfer 
measurement considered in this paper is for the same radiused 
edge blade in [3] but is equipped with a mean-camberline strip 
as given in Bunker and Bailey [4]. 

Figure 1 shows the definition of the airfoils and the shroud. 
The blade profiles are typical of a large power generation 



turbine. The cascade is linear and the span is 10.16cm. As can be 
seen the shroud incorporates a recess ahead of the blades to model 
a similar feature found in an actual turbine shroud. The 
instrumented blade has a rounded edge of 2.54-mm radius. The tip 
clearance addressed in this paper is the nominal value of 2.03 mm 
measured from the top surface of the mean-camberline strip. The 
height of the mean-camberline strip was 1.27 mm. This places the 
blade tip surface at 3.30 mm from the shroud (In [3] this distance 
was 2.03 mm). The strip in the experiment and in the calculations 
has a square cross section and sharp edges. Figure 2 shows the 
design of the actual two passage blade cascade as reproduced 
from Bunker et al. [3]. The extent of the strip is from the leading 
edge to 85% axial chord. Table 1 lists the run conditions for the 
cascade and input to the numerical simulations. 

THE COMPUTATIONAL METHOD 

The simulations in this study were performed using a multi- 
block computer code called Glenn-HT [6] which is based on a 
single block code designed by Amone et al. [7]. This code is a 
general purpose flow solver designed for simulations of flows in 
complicated geometries. The code solves the full compressible, 
Reynolds-averaged Navier-Stokes equations using a multi-stage 
Runge-Kutta based multigrid method. It uses the finite volume 
method to discretize the equations. The code uses central 
differencing together with artificial dissipation to discretize the 
convective terms. The overall accuracy of the code is second 
order. The present version of the code [5,8,9] employs the k-co 
turbulence model developed by Wilcox[10,ll] with modifications 


TABLE 1: Run Conditions 


Pressure ratio across the blade row 

0.69 

Exit Reynolds number 

2.57E6 

Inlet Mach number 

0.30 

Turbulence intensity 

5% 

Inlet angle 

44.9 deg. 
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by Menter [!2|. The model integrates to the walls and no wall 
functions are used. For heat transfer a constant value of 0.9 for 
turbulent Prandtl number, Pr^ and a value for Prandtl number 
(Pr) equal to 0.72 is used. Viscosity is a function of temperature 
through a 0.7 power law [ 1 3 land Cp is taken to be a constant. 

GEOMETRY MODELING AND GRID SYSTEM 

As in |21 the complete two blade passage was modeled. 
This was done to achieve better agreement with the measured 
pressure distribution on the tip in the cited work. Figure 3(a) 
shows the blade surface and the modeled mean-camberline strip. 
Figure 3(b) (reproduced from 12]) shows the grid for the 
complete passage including the splitter plate. Figure 3(c) shows 
the surface grid on the blade tip with the radiused edge and the 
mean- camberline strip. The grid details on the tip including the 
mean-camberline strip is shown in figure 3(d). 3(e) shows 
further the details of the grid around the strip. The topology is 
designed such that the viscous gridlines stay near the surfaces 
(including the strip) and do not spread in the “inviscid’' regions 


of the flow. The main differences between the present simulation 
and the actual geometry are twofold: The mean-camberline in the 
experiment is made up of piecewise straight lines and starts right 
at the leading edge whereas in the CFD simulation the strip is 
placed as a continuous strip and starts at a location 5^/r of axial 
chord and not from the leading edge. This was done for practical 
grid generation reasons. The grid is generated using a 
commercially available computer program called GridPro™. Tbe 
model of Fig. 3 consisted of 1 .2 million cells when a a sharp edge 
tip was u.sed. The number of cells were 1.4 million for the 
geometry with the radiused edge smooth tip. For the present case 
of a blade with a mean-camberline strip 1.8 million cells were 
used. The viscous grid is generated by embedding grid lines 
where needed, including the grid around the splitter. The 
stretching ratio did not exceed 1.25 for the viscous grid away 
from the no-slip surfaces. The distance to the first cell center 
adjacent to solid wall is such that the distance in wall units, (y"^) is 
near or below unity. Specifically in the tip region the number of 
grid cells from the tip of the strip to the shroud was 50 and from 
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Fig. 4 Streamline patterns of the flow in the tip region for 
(a) radius edged tip (b) radius edged tip and meaneamber 
line strip 


the tip surface to the shroud was 58 cells. Further details about 
the grid may he found in [2]. 

RESULTS AND DISCUSSION 

Along with the conditions given in Table 1, a very thin 
boundary layer thickness of 0.1% of the passage height was 
imposed at the inlet to the computational domain upstream of the 
flow splitter (at the location of upstream turbulence grid in Fig. 
2). The tunnel side walls were modeled using a slip boundary 
condition [2]. 

Figure 4 shows the streamline patterns of the flow in the lip 
clearance. Due to the rounded edge of the tip, the flow has little 
difficulty negotiating the turn into the clearance passage. In 
contrast on a sharp edge tip, a separation vortex forms along the 
pressure side edge of the blade tip thus narrowing the passage 
and thereby reducing the leakage[2]. Figure 4(b) shows the 
streamline pattern with the strip present. A separation vortex 
pattern behind the strip can be discerned from that figure. Figure 
5 shows the pressure distribution over the blade tip. Note that the 
areas of low pressure are associated with high velocities. In Fig. 
6, the velocity vectors over the blade tip for several axial cross 
sections (Fig. 6(a)) and a plane parallel to lip at the half height of 



Fig. 6 Velocity vectors colored by the velocity magnitude 

|v|. 
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Tip treatment 

Percent 
mass flow 

Y 

Radiused edge 

2.5% 

0.118 

Radiused edge+MCL Strip 

^ 2.0% 

0.125 

Sharp edge 

1.8% 

0.114 


TABLE 2: Tip leakage rates and total pressure losses 


the strip are plotted (Fig. 6(b)). The figures show that the flow 
behind the strip possesses high velocity component parallel to 
the strip. It also shows that the velocity is reduced as the strip is 
approached. This flow pattern has implications for the blade tip 
heat transfer which will be discussed later in this paper. 

Tip leakage rate and efficiency comparisons 

Tip leakage flow rate is often a good indicator of tip losses 
and tip heat transfer in relation to different tip treatments. 
Usually the smaller the flow rate the lower the losses and smaller 
the tip heat transfer. Since it is difficult to measure the tip 
leakage flow rate, the computations can serve to provide such 
information as has been done here. Table 2 below shows the tip 
leakage as a fraction of the total mass flow into the blade 
passage. Note that the total mass flow is for two passages. 

Table 2 shows that the mean-camberhne strip reduces the tip 
leakage by 20% as compared to the radiused edge blade tip 
without the strip. Interestingly for a sharp edge tip with the same 
clearance, as shown in that table, the reduction in the leakage 
mass flow rate is 28%. Also shown in Table 2 is the Y factor 
which is a measure of total pressure loss. The case with the 
mean-camberline strip although with lower tip leakage than the 
rounded case produces a larger loss. This also was the case with 
the squealer tip case of reference [5] where the reduction in 
leakage was not accompanied by a calculated improvement in 
efficiency. 

Heat Transfer 


The rate of heat transfer is presented in terms of heat 
transfer coefficient which is defined as: 


^wall 

wall ~^inlet^ 


Tinlet is the inlet total temperature and Twall is the wall 
temperature. Twall was determined from the experimental 
, measurements to be 1.06 times the inlet total temperature. The 
heat transfer coefficient was calculated on the blade tip where a 
constant temperature boundary condition is imposed. An 
effectively adiabatic boundary condition was imposed on all the 
other surfaces. 


In Fig. 7, the radiused edge blade tip heat transfer is 
reproduced from [2]. The agreement between the experiment 



Fig. 7 Simple radiused edge blade tip. Heat transfer 
coefficient for 2.03 mm clearance and Tu=5%. (W/m2/K), 
(a) Measured (Bunker et al, 2000) and (b) calculated 



Fig. 8 Radiused edge blade tip with mean-camberline strip. Heat 
transfer coefficient for 2.03 mm clearance and Tu=5%. (W/m2/ 
K), (a) Measured (Bunker and Bailey, 2000) and (b) calculated 


NASA/CR— 2001-210764 


5 




and simulation for this case is consistently better than 15% over 
the entire tip. Figure 8 shows experimental and calculated tip 
heat transfer for the tip with the mean-camberline strip. The area 
of the largest heat transfer coefficient appears to be midchord 
near the suction side on the exp>eriment and the simulation. 
Referring to Fig. 5, this area corresponds to the area of low 
surface pressure and thus high velocity over the blade tip. There 
are some discrepancies between the simulation and the 
experiment. They include the magnitude of the low island on the 
suction side near the strip and the apparent low value adjacent to 
the strip on the pressure side. These are likely caused by 
geometry differences between the model and the experiment 
discussed earlier or by the inadequacies of the turbulence model. 
The breakdown of the one-dimensional conduction assumption 
for the experimental data near the edges can also play a role in 
those differences. In the simulations the rate of heat transfer 
drops as the strip is approached on the pressure side of the strip 
and decreases away from the strip on the suction side. Figure 6 
shows that the fluid velocity is diminished as it approaches the 
strip on the pressure side, this reduction in velocity reduces the 
heat convection to the blade. On the suction side of the strip there 
is a vortical zone within which the flow has a large velocity 
component parallel to the strip. This streamwise vortex 
continually feeds core temperature air to the blade tip, thus 
increasing the heat transfer. This accounts for the high heat 
transfer rate predicted by the numerical simulation. 

It should be noted that without the mean-camberline strip, 
due to the increase in the clearance, the rate of heat transfer for 
the case of Fig. 8 would be larger than that of Fig. 7 for most of 
the blade tip surface except near the tail. The tail area in Fig. 6 
does have a smaller rate of heat transfer. This has been reported 
in [14]. 

SUMMARY AND CONCLUSIONS 

In this paper the numerical prediction of flow and heat 
transfer and comparison with the experimental heat transfer data 
of tip heat transfer for a blade equipped with a tip mean- 
camberline strip (squealer) was perf^ormed. The blade chosen 
represents a first stage blade of a large power generation turbine. 
The casing upstream of the blade tip was recessed. The effect of 
mean-camberline strip on the tip leakage flow was quantified and 
compared to plain tip surfaces with radiused edge and sharp 
edges. It was found that the sharp edge works better than the 
radiused edge in reducing the tip leakage flow. For the radiused 
edge case it was found that the mean-camberline strip reduces 
the tip leakage but the total pressure loss is not commensuralely 
reduced. Reasonably good comparison with the experimental 
heat transfer was achieved. Areas of large heat transfer rate were 
identified and likely reasons for local heat transfer enhancement 
was given using flow patterns. 
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